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ABSTRACT 

We conducted radio interferometric observations of six pulsars at 610 MHz using the Giant 
Metrewave Radio Telescope (GMRT). All these objects were claimed or suspected to be the 
gigahertz-peaked spectra (GPS) pulsars. For a half of the sources in our sample the interfer¬ 
ometric imaging provides the only means to estimate their flux at 610 MHz due to a strong 
pulse scatter-broadening. In our case, these pulsars have very high dispersion measure values 
and we present their spectra containing for the first time low-frequency measurements. The 
remaining three pulsars were observed at low frequencies using the conventional pulsar flux 
measurement method. The interferometric imaging technique allowed us to re-examine their 
fluxes at 610 MHz. We were able to confirm the GPS feature in the PSR B1823 —13 spectrum 
and select a GPS candidate pulsar. These results clearly demonstrate that the interferometric 
imaging technique can be successfully applied to estimate flux density of pulsars even in the 
presence of strong scattering. 

Key words: pulsars: general - pulsars: individual: B1750—24, B1800—21, B1815—14, 
B1822-14, B1823-13, B1849H-00 


1 INTRODUCTION 

In the case of most pulsars, their observed radio spectra can be de¬ 
scribed using a power law with a negative spectral index of —1.8 
or (for a small fraction of sources) two power laws with spec¬ 
tral indices of —0 .9 and —2.2 with a break frequency i/b on av¬ 
erage of 1.5 GHz iMaron et~n]|2000h . Some pulsars also exhibit a 
low-f requency turnover in their spectra jSiebeJl973l : lLorimer et al.l 
Il995t) . A spectrum of that kind is characterized by a positive spec¬ 
tral index below a peak frequency Vp of about 100 MHz (with a 
few exceptions when the spect rum peaks at frequen cies up to sev¬ 
eral hundred MHz). However, iKiiak et ^ bOl Ibl) pointed out a 
small sample of pulsars that peak around 1 GHz and above. Such 
an object, called the gigahertz-peaked spectrum (GPS) pulsar, is de¬ 
scribed as a relatively young source that has a high dispersion mea¬ 
sure (DM) and usually adjoins a dense, sometimes extreme vicinity. 
This suggests that the GPS in pulsars might be caused by either the 
conditions around neutron stars or the physical properties of the 
interstellar medium. 

The strongest argument for environmental origin of the high- 
frequency turnover in radi o pulsars spectra is the evolution of 
PSR B1259—63 spectrum. iKiiak et ^ j2011ah showed that the 
spectrum of the pulsar at the various orbital phases exhibits both 
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a shape and a peak frequency evolution due to the orbital motion 
of the pulsar around its companion Be star LS 2883 on a very el¬ 
liptical orbit. The PSR B 1259—63 spectrum demonstrates a strong 
similarity with the gigahertz-peaked spectra, especiall y when the 

S iulsar i n its motion gets closer to its companion star. IKiiak et aH 
201 T3) proposed two effects which can be responsible for the ob¬ 
served variations, the free-free absorption in the stellar wind and 
the cyclotron resonance in the magnetic field associated with the 
disk of Be star. Both these processes assume the absorption to be 
caused by external factors , like in the cases of the isolated GPS 
pulsars jKiiaketal .11201 lal) . 

IKiiak et H] l l2013l) studied the radio spectra of two magnetars 
PSRs J1550—5418 and J1622—4950 and in both cases their radio 
spectra clearly peak at the frequencies of a few GHz. Both these 
magnetars are associated with supernova remnants and hence sur¬ 
rounded by ionized gas which can be responsible for the free-free 
absorption of the radio waves. The authors concluded that the GPS 
feature in radio magnetars spectra can be of environmental origin, 
in the same way as it occurs in the vicinity of GPS pulsars. 

Pulsars with a high-frequency turnover in their spectra 
have represented th e smallest group of the r a dio p ulsar spec¬ 
tra types. However, [Bates. Lorimer & VerbiestI j2013l) estimated 
that the number of such sources may constitute up to 10% of 
the whole pulsar population. The sample of GPS pulsars was 
extended to include PSR J2007-I-2722, whose flux density mea- 
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sureme nts were presen t ed bvlAllen. Knispel & Corded ( l2013h . Re¬ 
cently, IPembska et alj ( l2014h reported two newly-identified GPS 
pulsars. One of them, PSR B1740-1-1000, is the first low-DM pulsar 
that exhibits the gigahertz-peaked spectrum. This case, along with 
the GPS phenomenon in radio magnetars, led the authors to con¬ 
clude that the GPS candidate selection criteria need to be revisited. 
In future searches for new GPS pulsars, the presence of interesting 
(or extreme) environments, instead of the high DM, could play a 
crucial role in the source selection process. 

IPembska et alj ( l2014h also pointed out that the small num¬ 
ber of the currently known GPS pulsars may be the result of our 
limited knowledge of pulsar spectra in general, especially below 
1 GHz. The authors outlined the need for a more extensive sam¬ 
ple of GPS sources to establish a plausible statistics about those 
objects. However, in the cases of some GPS candidate pulsars the 
standard pulsar flux measurement methods are affected by strong 
scattering at low frequencies. The phenomenon causes the pulse 
profiles to become broader, i.e. pulses attain roughly exponen¬ 
tially decaying scattering tail. It has been shown that the charac¬ 
teristic broadening of the pulse, Tsc, depends on both the observ¬ 
ing freguenc};j_as_welI_asDM_^e empirical relation was given 
by [Bhat. Cordes & C amillol 2004b. Recent results on scattering 
were discussed bv lLewandowski et aP j2013l) in their analysis of 45 
pulsars, based on the Giant Metrewave Radio Telescope (GMRT) 
and the Effelsberg Radio Telescope observations. Since the scatter¬ 
ing becomes stronger at lower frequencies, for a given pulsar the 
flux becomes increasingly underestimated then. For high-DM pul¬ 
sars at low frequencies, when the scattering time is greater than the 
pulsar period by a significant factor, one will see no pulsed emis¬ 
sion. Thus, the flux density measurements required to construct ra¬ 
dio pulsar spectra using “traditional” methods can be difficult or 
sometimes impossible to conduct. For these cases the only way to 
determine the pulsar flu x is using the interfer ometric imaging tech¬ 
niques (see for example lKouwenhovenll2()00h . 

The interferometric measurements of pulsar fluxes at both 
325 MHz and 610 MHz usin g the ( 3MRT have been demonstrated 
in iBasu. Athreva & Mitral OOllh and [Basu. Mitra & Athreval 
(l2012h . The imaging techniques provide a superior alternative to 
the standard flux measurements, especially in our studies since 
the sources we selected for observations are high-DM pulsars. For 
some of them the imaging techniques are the only secure means 
to estimate their flux. There is at least two reasons for employing 
imaging techniques. Firstly, flux calibration in an interferometer is 
more robust due to the baseline lying at zero level thereby reduc¬ 
ing errors made during the baseline subtraction of a normal pul¬ 
sar observation. Secondly, the instrumental and atmospheric gain 
fluctuations on very short time scales can be corrected using self¬ 
calibration of the interferometric data. The corrections are deter¬ 
mined by flux densities of constant and bright background sources 
in the field and hence would not be affected by the pulse variation 
of the relatively weak pulsar at a field center. 

In this paper we present flux measurements of six pulsars ob¬ 
served at the 610 MHz frequency band of the GMRT using inter¬ 
ferometric imaging technique. The sources selected for our studies 
are GPS pulsars or GPS candidates. We chose the 610 MHz band 
for these studies, firstly, because the frequency is low enough to es¬ 
timate whether a given object is indeed a GPS pulsar, secondly, due 
to the higher probability of detecting GPS pulsars at 610 MHz than 
at lower frequencies as a result of their inverted spectra at sub-GHz 
frequencies and finally, to avoid the RFI and other systematic ef¬ 
fects that are prominent at lower frequencies. This analysis allowed 
us to inspect flux densities of some objects and select a strong GPS 


Table 1. The list of sources and their parameters, where DM - dispersion 
measure, P is period and Tsc denotes scattering time. 


Pulsar 

DM 

(^) 

Age 

(kyr) 

P 

(sec) 

Tsc 

(sec) 

Associations 

B1750-24 

672 

593 

0.528 

2.38 


B1815-14 

622 

2270 

0.291 

0.457 


B1849+00 

787 

356 

2.18 

6.19 


B1800-21 

234 

15.8 

0.133 

0.039 

1,2,3 

B1822-14 

357 

195 

0.143 

0.279 

1 

B1823-13 

231 

21.4 

0.101 

0.0337 

1,2 


1 - HESS (High Energy Stereoscopic System) 

2 - X-PWN (X-ray pulsar wind nebula), 3 - SNR (supernova remnant) 
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Figure 1. The figure shows the variation of the flux value of the surrounding 
sources in the field of view across the three obseiving runs. The strong 
sources in the fields with the pulsars PSRs B1750—24 and B1849+00 at the 
field centers were identified and their average flux across the three sessions 
were calculated. The ratio of the flux of each source with respect to the 
average flux and as a function of the average flux is plotted in the figure, 
with each session marked with different symbols. The ratios are scattered 
around unity with the noise in the scatter decreasing with the increasing flux 
levels. This demonstrates the analysis process to be correct and consistent 
for all the three observing sessions. 


candidate pulsar. We were able to confirm that the the interferomet¬ 
ric imaging technique can be successfully applied to estimate flux 
density of pulsars. 


2 OBSERVATIONS AND DATA ANALYSIS 

We recorded the interferometric data on six pulsars using the 
GMRT that is located near Pune, India and consists of an array of 
30 distinct dishes, each with a diameter of 45 meters, and a total of 
435 baselines. The dishes are spread out over a region of ~ 27 km 
and roughly resemble a Y-shaped array. The data were recorded at 
the 610 MHz frequency band with a bandwidth of 33 MHz spread 
over 256 frequency channels. All the six sources were observed on 
three separate dates: 30 December 2012, 6 January 2013 and 13 
January 2013, each observation separated by a week, to account for 
variations in the pulsar flux over long time scales. 

The observations were carried out using standard schemes 
where strategically placed calibrators were interspersed with the 
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Table 2. Flux density measurements resulted from the interferometric imag¬ 
ing observations in three observing sessions (5i, S 2 and S 3 respectively) 
along with the uncertainties which include calibration errors, the rms noise 
in the maps and fitting en'ors. (5) denotes the weighted mean of all results 
(presented along with its uncertainty) for a given pulsar. 


Pulsar 

Si 

(mJy) 

S2 

(mJy) 

53 

(mJy) 

(S) 

(mJy) 

B1750-24 

3.65T0.28 

3.90T0.25 

3.98T0.25 

3.9T0.3 

B1815-14 

24.5T1.8 

24.6T1.5 

26.1T1.5 

25.1T1.6 

B1849-1-00 

15.2T1.1 

15.5T1.0 

16.1T1.0 

15.39T0.95 

B1800-21 

8.22T0.60 

6.75T0.42 

7.82T0.48 

7.4T0.9 

B1822-14 

3.31T0.36 

3.43T0.35 

4.18T0.38 

3.5T0.6 

B1823-13 

3.51T0.26 

3.42T0.22 

3.53T0.23 

3.5T0.2 


sources. We recorded flux calibrators 3C48 and 3C286 before and 
after each observing run for 10 min. The antenna gain variations 
were determined using a phase calibrator observed every half hour. 
The sources in our sample are nearby in terms of their position in 
the sky and the phase calibrator 1822—096 was used for all the 
sources to carry out amplitude and phase calibrations as well as 
bandpass calibration to correct for variation across frequency band. 
Each of the pulsars was observed for 10 minutes at a time and six 
such observations on each source were well spaced out during each 
of the observing sessions which lasted for about eight hours. We 
recorded one hour of data on each pulsar during each observing 
run with all the sources having a wide coverage in the u-v plane 
due to our strategy of spacing out their pointings. 

The data reduction was done using the Astronomical Image 
Processing System (AIPS) where we used standard techniques of 
RFI flagging, calibration and imagi ng. The flux scales were deter¬ 
mined using the measurements of [Baars et al.l l ll977h with latest 
corrections by Perley et al. 2010 (in AIPS). All the pulsar fields 
during each observing run was passed through multiple rounds of 
self-calibration and imaging to improve noise characteristics in the 
maps. The following checks were performed to test the fidelity of 
the images: 

1. The flux of the source 1822—096 (phase calibrator) was deter¬ 
mined for each of the three observing sessions using both the flux 
calibrators 3C48 and 3C286. The flux value for each of the sessions 
were 6.3±0.5 on 30 December 2012, 6.2±0.5 on 6 January 2013 
and 6.3±0.5 on 13 January 2013. This shows that the various inde¬ 
pendent calibrations applied during each separate observations and 
also different flux calibrators gave identical results demonstrating 
the technique to be correct. In addition the correctness of the flux 
scale was demonstrated by comparing the flux value of the phase 
calibrator to the expected value (between 325 MHz and 1.4 GHz) 
which yielded consistent results. 

2. We performed an additional test for determining the correctness 
of the analysis techniques. The same fields were observed on mul¬ 
tiple days with independent calibrations for the system. Each of 
the pulsar fields imaged, with the pulsar at a field center, had a 
half power radius of ~ 40 arcmin with a large number of sources 
(> 100) in each field of view. If an analysis is correct then the flux 
of all these sources, in an ideal scenario, should be identical for all 
the observations. However, under realistic conditions one expects 
some noise across the different observing runs. To demonstrate this 
we used two (out of six) pulsar fields and determined the flux val¬ 
ues of all the strong sources (flux greater than 10 times the rms 
noise in a map) for all the three observing runs. An average flux 


for each source across the three observing runs was calculated and 
the ratio of the flux to average flux was plotted as a function of the 
source flux value (see Figure[TJ. As is clear from the figure the flux 
ratios are concentrated around unity with the noise increasing with 
decreasing flux level of a source. This is exactly as per our expec¬ 
tations and clearly demonstrates the correctness of our calibration 
and analysis techniques across the three observing sessions. 


3 RESULTS 


Table [T] gives the list of sources with some of their basic pa¬ 
rameters (DM, age and period). All these objects are either con¬ 
firmed or candidate GPS sources. As it is clear from the table, the 
sample was subdivided into two categories. We have made rough 
estim ates of the scattering tim escales using either observational 
data jLewandowski et ani2013h or the predictions derived from a 
single thin screen model (using the scatter time frequency scaling 
index a = 4). For the first group of pulsars, PSRs B1750—24, 
B1815 —14 and B1849-I-00, the pulse broadening at 610 MHz is 
large, hence the interferometric imaging provides the only mean 
to estimate flux for these sources. The flux density of the re¬ 
maining pulsars, PSRs B1800—21, B1822—14 and B1823 —13, 
was me asured at low frequencies by conventional pulsar obser¬ 
vations I Lorimer et alj 1 19951 : iKiiak, Gupta & Krzeszowskil 1 20071 : 
iKiiak et alj201 Ibl) . however since the predicted (or measured) scat¬ 
ter time estimates are a significant fraction of the pulsar periods, 
there is a possibility for the flux densities measured that way to be 
underestimated. Additionally, these three sources have counterparts 
in the High Energy Stereoscopic System (HESS) obser vations with 
indicati ons o f pulsar win d nebula (PWN) around them. IKiiak et al.l 
(l2011ah and IKiiak et Zl ( 1201 Ibl) reported PSRs B1822-14 and 
B1823 —13 as GPS pulsars, and PSR B1800—21 is considered as 
a GPS candidate. We included these pulsars in our studies to re¬ 
examine their flux values at 610 MHz where some profile broaden¬ 
ing due to scattering is also present. 

For the six pulsars we constructed radi o spectra, combinin: 
flux density measurements from the literature |L orimCTetalJ_199; 


Maron et^l200d : Kiiak, Gupta & Krzeszowskill2007l : IKiiak et H] 


201 Ibl . and references therein) and the ATNF (Australian Tele¬ 


scope National Facility) pulsar catalogu^ together with the new 
results shown in Table (see online material for the maps of 
the pulsars that were used to estimate their flux density at 610 
MHz). The spectra are presented in groups, depending on their 
morphological properties. Our studies show that for the high-DM 
pulsars, namely PSRs B1750—24, B1815—14 and B1849-I-00, 
their spectra resemble a simple power law. We were able to 
confirm the GPS feature in the spectrum of PSR B1823—13 
and pointed out two sources, PSRs B1822—14 and B1800—21, 
suspected to be gigahertz-peaked spectra pulsars which require 
further investigation. The fits presented in the paper were obtained 
by the implementation of the nonlinear least-squares Levenberg- 
Marquardt algorithm.The results of the fitting procedure are given 
in Tab.[3 


PSRs B1750-24, B1815-14 and B1849-t00 have very high 
DMs in excess of 600 pc cm~®. Their spectra are shown in Figurej^ 
These objects are characterized by a significant pulse broadening 
due to interstellar scattering (of the order of several periods). As 


http://www.atnf.csiro.au/research/pulsar/psrcat/ iManchester et alj2005h 
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Figure 2. The spectra of PSRs B1750—24, B1815—14 and B1849+00, pul¬ 
sars with high DMs (see Tab. \Q. Open circles denote the GMRT inter¬ 
ferometric obser vations. Measureme n ts marked with black dots are taken 
from liter ature jMaron et a 0 l2QQ(t iKiiak. Gupta & Krzeszows'O l2007l : 
iKiiak et ^12011 J and references therein) and the ATNF pulsar catalogue. 
The straight lines represent our fits to the data using power-law function 
(in the case of PSR B1815—14 the flux measurements at 610 MHz marked 
with black dots were excluded from the fitting procedure). The power in¬ 
dices resulted from the fitting procedure are given on each panel, for more 
details see Tab.[^ 


mentioned above, the sources were suspected to be GPS pulsars 
but our analysis suggest that their spectra can be fit by a power law. 

In the case of PSR B1815—14 the interferometric imaging 
technique was used to inspect its flux density measurements at 610 
MHz due t o a strong pulse scatter-b roadening (over a full pulse 
period, see iLewandowski et al.ll201^ . Previous observations sug¬ 
gested that the pulsar spectrum shows a turnover feature. Our new 
result, accompanied by a large scatter time estimate, suggests a sig¬ 
nificant underestimation in the previous flux measurements for this 
pulsar at 610 MHz. 


Table 3. Fitted parameters to data of four pulsars, PSRs B1750—24, 
B1815—14 and B1849+00 (using a power-law, where g is a power 
index) and B1823—13 (using the same function (1) as in IKiiak et al.l 
(201 lah ). The implementation of the nonlinear least-squares Levenberg- 
Marquardt algorithm was used to perform the fitting procedure. The pa¬ 
rameters are given with a reduced ■ For more details on data excluded 
from the fitting procedure see Figs[2and[^ 


Pulsar 

Fitted parameters 


B1750-24 

€ = -1.0 ±0.14 

7.5 

B1815-14 

€ = -1.75 ±0.07 

0.95 

B1822-14 

^ = -0.48 ±0.08 

2.5 

B1849+00 

€ = -1.9 ±0.2 

9.6 

B1823-13 

a = -0.95 ± 0.18, b = 0.28 ± 0.11 
c = 0.61 ±0.03 

3.4 


The classification of the PSRs B1750—24 and B1849+00 
spectra was not possible before due to very limited flux density 
measurements at frequencies below 1 GHz. It seems clear that the 
spectrum of PSR B1849+00 one is a typical steep pulsar spectrum. 
The PSR B1750—24 spectrum can also be described by a single 
power law but at the same time it is relatively flatter than a usual 
pulsar spectrum (a spectral index of —1.0; the fits and the resulting 
spectral indices are presented in Fig.^. 

PSR B1823—13, whose spectru m is presented in Fig ure 
was classified as a GPS pulsar by iKii^ et ^ (201 Itj) . The 
object was included in our sample to re-examine its flux density 
measurement at 610 MHz. It is noteworthy that new measurements 
indicate a larger flux density than previous GMRT observations 
which were carried out using the instrument in phased array 
mode, t hus standard pulsa r flux measurement methods were 
applied (Kiiak et alj 120119) . Even if one disregards the earlier 
measurements and uses only the interferometrically derived flux, 
the spectrum continues to exhibit a GPS feature. However, it seems 
clear that this object requires further observations at frequencies 
below 600 MHz which we plan to perform in future projects. 

PSR B1800—21, whose spectrum is presented in Fig. |4l can 
be treated as a new, very promising GPS candidate pulsar. It 
is a young, Vela-like p ulsar, associated with a supernova rem¬ 
nant dKiiak et alj201 IbO . Previously the spectrum seemed to be flat 
at low frequencies but our estimate, which indicates a much smaller 
flux density than in the previous measurements at 610 MHz, sug¬ 
gests a positive spectral index in the low-frequency range. New re¬ 
sults clearly imply that the earlier flux measurements should be ver¬ 
ified. 

PSR B1822-14 is the last pulsar in our sample. Its spectrum 
is presented in Figure|5] Similar to the case of PSR B1823—13, the 
flux density measured using the interferometric imaging techniques 
is greater than the values from s tandard method. PS R B1822—14 
was identified as a GPS pulsar (Kiiak et al.ll201 Ibh . however the 
interferometric measurements suggest otherwise - the spectrum of 
the pulsar seems to be a power law (with spectral index of —0.48) 
when including new measurements which may indicate that the 
earlier measurements were affected by interstellar scattering. This 
pulsar definitely needs further investigation and possibly additional 
measurements at lower frequencies. 
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PSRB1823-13 (J1826-1334) 


PSRB1800-21 (J1803-2137) 




Figure 3. The spectrum of the GPS pulsar PSR B1823—13. Open circle 
denotes the GMRT i nterferometric observations, whereas black dots denote 
data from literatur e iMaron et alj2003:lKiiak. Gupta & Krzeszowskill2007l : 
iKiiak et'ani2Ql laL and references therein) and the ATNF pulsar catalogue. 
T he curve represents our fits to the data using the same function (1) as 
in iKiiak et^j201 ls^ ■ The measurement at 610 MHz marked with the black 
dot was excluded from the fitting procedure. For fitted parameters see Tab.|3] 


4 DISCUSSION 

The interferometric imaging technique allowed us to measure flux 
density for six pulsars with the GMRT at 610 MHz. The frequency 
was chosen to verify whether the objects are GPS pulsars. In our 
studies we included three pulsars which had previous flux measure¬ 
ments at 610 MHz (obtained by using the standard techniques) in 
order to re-examine these estimates. 

In the case of the remaining three sources, due to a strong 
pulse broadening, the interferometric observations are the only way 
to estimate their flux at this frequency (see the scatter time esti¬ 
mates in Tab. [T). We inspected the spectra of these three pulsars 
with very high-DMs (>600 pc cm"®). PSRs 1750-24, 1815-14 
and B1849-1-00 were suspected to be GPS sources but with the addi¬ 
tion of new measurements their spectra appear to resemble a typical 
single power-law pulsar spectrum. We can also confirm the GPS 
feature in the PSR B1823 —13 spectrum and select a new strong 
GPS candidate: PSR 1800—21. Its spectrum as well as the spectrum 
of PSRB1822—14, require however some further investigation. 

4.1 The interstellar scintillations and scattering 

The variations in flux density values resulted when applying dif¬ 
ferent observing techniques can be caused by interstellar scintilla¬ 
tions. The pulsars in our sample have DMs ranging from ~200 to 
--^800 pc cm"® which means that any scintillations will occur in a 
strong scintillation regime. Hence we have to take into considera¬ 
tion both diffractive and refractive scintillations. As the diffractive 
time scales, Atmss oc where d denotes distance, in a 

case of high-DM and hence distant sources should be up to a few 
minutes at the observing frequency / of 610 MHz, thereby a single 
1 hour observing run should be enough to averaged-out any varia¬ 
tions of the measured flux density. The discrepancies with the ear¬ 
lier flux density measurements (if present) can be at least partially 
caused by the refractive interstellar scintillations, as the refractive 
time scales, Afmss oc /"^'^d^ ®, may vary from days to months. 
However, the refractive modulation index Ammss oc jO. 57 ^- 0.37 


Figure 4. The spectrum of the candidate for GPS pulsar, PSR B1800—21. 
Open circle denotes the GMRT intei-ferometric obse rvations. Measure¬ 
ments marked with black dots are taken from literature jMaron et alibood : 
IKiiak et alJ201 1 j. and references therein) and the ATNF pulsar catalogue. 


for high-DM pulsars should be relatively small at the frequency 
of 610 MHz. To account for refractive scintillations influence, we 
performed our flux measurements on three epochs (separated by no 
less than a week). For the purpose of constructing pulsar spectra we 
are using an average of values obtained from these measurements. 
We have to note however, that in this regard the interferometrically 
measured flux density values are not affected any more than the 
ones resulting from the standard observing procedure. Pulsar flux 
measurements regardless of the method used will be affected by the 
scintillation-driven flux variation in exactly the same way. 

As mentioned in the Introduction, for high-DM pulsars, and 
especially when the observations are conducted at low observing 
frequencies, the pulsar profiles can be affected by interstellar scat¬ 
tering to such an extent that it is difficult, or even impossible to per¬ 
form flux measurement using the standard profile-based method. 
Strong scattering causes the profiles to attain scattering tails, what 
may significantly affect the observed profile background level so 
that the proper baseline required for the standard flux measurement 
can not be found. Especially prone to such errors will be the cases 
where the scatter time is comparable to the pulse period, as one 
will still be able to see a prominent pulse, while the scattering tail 
will hide the proper baseline level. One has to remember that even 
a “moderate” length of the scattering tail may still affect the profile 
baseline and cause erroneous results, es pecially when the observa - 
tions are affected by a high noise level JPewandowski et ^l2013h . 
In some circumstances the influence of scattering may even change 
the appearance of the spectrum by mimicking a turnover feature, 
which may lead to an erroneous identification of a pulsar as a GPS 
source, like in the case of PSR B1815 —14. 

Since there is no simple way to account for the scattering- 
induced baseline level change one should refrain from using the 
standard method in such cases. Practically, however, it can be hard 
to judge if we are dealing with a profile where scattering led to an 
erroneous flux density estimation. We believe that at least some of 
the past flux measurements performed for the pulsars from our sam¬ 
ple might have been affected by the scattering, which is discussed 
below. Regardless, to avoid the scattering-induced issues one needs 
a reliable method of pulsar flux density measurement that will not 
be affected by scattering. 

















6 M. Dembska et al. 


4.2 Imaging technique and pulsar flux density 

Interferometric imaging techniques at low radio frequencies have 
been used to determine flux density of a wide class of astronomical 
sources and can be easily extended to determine pulsar flux values. 
The estimated flux of the phase calibrator 1822—096 at 610 MHz 
was identical over the three widely separated observing sessions 
and using two different flux calibrators 3C48 and 3C286. In ad¬ 
dition its flux estimates are in the expected range when compared 
with the flux densities at 325 MHz and 1.4 GHz (from the Very 
Large Array calibrator manual). The consistency of our calibrations 
over the different observing sessions was also demonstrated by the 
flux ratios of the surrounding background sources in the field of 
view, which hovered around unity (see Figure[T](. Additional vindi¬ 
cation of our results is provided by the flux density measurements 
of the three pulsars B1750—24, B1815 —14 and B1849-I-00 which 
could not be determined at 610 MHz using conventional pulsar flux 
measurement techniques due to their highly scattered profiles. Our 
flux estimates at 610 MHz for all the three cases are consistent with 
the pulsar spectra determined from higher frequencies, where inter¬ 
stellar scattering should not affect the measurement (see Figurej^. 

We now examine the cases of the three other pulsars 
B1800—21, B1822—14 and B1823 —13 where traditional pulsar 
flux measurements have been carried out in the past. All of these 
sources adjoin interesting environments which m ay cause an ad¬ 
dition al absorption at frequencies below 1 GHz. IPembska et aH 
( I 2 OI 4 I) suggested that an influence of such environment on pul¬ 
sar radio emission can manifest itself as a reverse spectrum with 
positive spectral index, which made these objects a plausible GPS 
candidates. 

Two of the pulsars, PSRs B1822—14 and B1823 —13, were 
classified by iKii^ et ^ ^201 lah and iKiiak et H] (1201 Ibt) as GPS 
pulsars. Obviously, one can note the discrepancies between stan¬ 
dard and interferometric flux measurements for both sources. The 
interferometrically measured flux densities are larger than the ones 
which resulted from standard observations. We know that the pulse 
pro file of PSR B1822—14 at 6 10 MHz shows a significant scatter¬ 
ing jLewandowski et ^l2013h . thus the flux density obtained us¬ 
ing the standard method is very likely to be underes timated. In the 
case o f PSR B1823 —13 the profile that was used bv iLorimer et al.l 
(Il995h to estimate the flux density of the pulsar at 610 MHz was 
not published. However, the estimated scatter time at this frequency 
is a significant fraction of the pulsar period (as in the case of PSR 
B1822—14). Hence one cannot rule out the possibility that conven¬ 
tional observations were in this case also affected by scattering. 


4.3 The case of PSR B1800-21 

As we pointed out above, usually the interferometric observations 
give better estimations of previously underestimated flux density 
values. Yet, one can note that in the case of PSR B1800—21 
its interferometrically measured flux value is lower that t he one 
res ulted from previous measurements by iLorimer et al. I 
and iKiiak et ^ ( 1201 Ibll . The discrepancy cannot be explained in 
a simple way and it is not possible to definitely point out the 
reason behind it. However, one still can propose possible effects 
causing such a variation. Assuming our latest result to be accurate 
and precise, there appears a chance that the difference is caused 
by rather external (i.e. environmental) factors. PSR B1800—21 is 
a young Vela-like pulsar, located at the southw estern edge of the 
G8.70 .1 nebula, a shell-type supernova remnant dKassim & Weiled 
Il990al f9). It has a plausible association with the 7 -ray source HESS 


PSR B1822-14 (J1825-1446) 



Figure 5. The spectrum of the GPS pulsar PSR B1822—14. Open circle 
denotes the GMRT i nterferometric ohseivations, whereas black dots denote 
data from literatur e iMaron et alj2QQQHKiiak, Gupta & Krzeszowsl3l2QQ7l : 
IKiiak et 511201 lal and references therein) and the AT NF pulsar ca t alogue. 
The dashed curve represents the fit to the data made bv iKiiak et alJi2011aL 
see their Fig. 1, a point marked with an open circle was not included). We 
present our power-law fit (the straight line) to the data, excluding two mea¬ 
surements at 610 MHz marked with black dots. The fitted power index is 
given with the fit, for more details on the fitting results see Tab. [3] 


J1804-216 dAharonian et al.ll 200 l | 200 ^: Iffigashi et alj|2008h . All 
these objects, together with a number of discrete ultracompact 
HII r egions, compose the W30 complex of peculiar morphol¬ 
ogy dFinlev & Oegelm^ Il994h . which emits over frequencies 
ranging from radio to 7 -rays. Thus, it is possible that during our 
interferometric observations the observed PSR B1800—21 flux was 
affected by an additional absorption, which manifested itself at 
lower frequencies. Such transient absorption may appear due to rel¬ 
ative motion of the components of the W30 complex. 

One should note that the two standard measurements which 
yielded higher flux values were separated in time (by more than 
10 years) and conducted independently using different observing 
facilities. Nevertheless, one cannot exclude the influence of the 
classical pulsar flux measurement procedure, such as non-linearity 
of a receiver or some non-typical flux calibration issues. Either 
way, PSR B1800—21 definitely requires further investigation. We 
have already conducted another set of interferometric observations 
of this source at 610 MHz with the GMRT to verify its flux density. 
The data is still to be analysed, however preliminary results show 
that the flux density resulted from our observations conducted 
on 28 December 2013 is 8.79 ± 0.67 mjy which is in a good 
agreement with values presented in Tab. Moreover, we have 
time allocation for 325 MHz observations of this source - these 
observations are to be conducted in January 2015. Consequently, 
we will be able to re-examine our results and hopefully resolve the 
discrepancy with the standard method measurements. 


5 SUMMARY 

The interferometric imaging technique allowed us to estimate 
pulsar flux density, making it possible to observe weak or high-DM 
pulsars at frequencies below 1 GHz and re-examine our previous 
results. In the cases of objects for which determining their flux at 
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low frequencies using standard pulsar observing methods is not 
possible (i.e. for a large Ts when comparing to pulsar period), 
it is the only way to estimate flux for these sources. Thus, the 
interferometric imaging technique provides a superior alternative 
to the pulsar flux measurements. We believe that the technique will 
help us to confirm more GPS pulsars in future. 

Moreover, we used this method to rule out three GPS pulsar 
candidates and confirm one GPS pulsar, PSR B1823 —13. We also 
pointed out two GPS candidates, including a very promising one, 
PSR 1800—14, that require further investigation. Hence, our analy¬ 
sis clearly shows that the interferometric imaging technique can be 
successfully applied to estimate flux density of pulsars. The method 
is the only secure way to determine the flux density of high-DM 
pulsars which are highly scattered at low radio frequencies. 
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